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Abstract 



Development is a complex phenomenon where the forming phenotype interacts with genetic 
and environmental inputs. Teeth are an important model for developmental studies and their 
development has been thoroughly investigated. However, despite of an extensive literature on 
the genetics of dental development, no studies have yet focused on the environmental influences 
on dental morphology. Here we aim to test whether and to what extent the environment plays 
a role in producing morphological variation in human teeth. We selected a sample of modern 
human skulls and used dental enamel hypoplasia as an environmental stress marker to identify two 
groups with different stress levels, referred to as SG ("stressed" group) and NSG ("non-stressed" 
group). We collected data on the occurrence and the relative development of 15 morphological 
traits (5 for each molar) on upper molars using a standard methodology commonly used in dental 
anthropology (ASU-DAS system) and then we compared the frequencies of the traits in the two 
groups. Overall, the results suggest that: (i) stressors like malnutrition and/or systemic diseases 
have a significant effect on upper molar morphology; (ii) stress generates a developmental response 
which increases the morphological variability of the SG; (iii) the increase in variability is directional, 
since individuals belonging to the SG have increased cusp dimensions and number. The se results 
are c onsistent with the expectations of the morphodynamic model of dental development ( Jernvalj 
1X9951) . 



Introduction 

In developmental studies, teeth represent an ex- 
cellent model, since once their development is 
completed, they undergo no further changes in 
size and shape (except from wear), thus record- 
ing the developmental process in their pheno- 
type. Our knowledge on the mechanisms of den- 



tal development has increased exponentially over 
the past 20 years. In particular, it has become 
evident how dental development, like the devel- 
opment of other organs, is regulated by inductive 



mal cells ( 


Ruch 




1995; 


Thesleff and Hurmerinta. 


1981 




Thesleff et al.. 


1995). The bulk of the 



research agenda has focused on t he genetic 
and molecular bases (reviewed in 
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Jernvall and Jung 20001 ): detailed knowledge is 
available today on over 300 genes expressed in 
mouse dental tissues and their patterns of ex- 
pression from the beginning of the tooth bud to 
the end of its development (see for example the 
Gene expression in tooth — WWW database, 
1996). The focus on the genetics of dental devel- 
opment has made it possible to make progress in 
our understanding of morphogenesis and varia- 
tion in dental morphologies; at the same time, 
it has overlooked other important aspects, such 
as the environmental influence on dental devel- 
opment and its contribution to the final pheno- 
type. It is probably because of this bias toward 
the genetic component of dental development re- 
search that is common to find in the literature 
such statements as "the devel opment of teeth 



is under strict genetic contro l " (iGalluccio et al 



2OI2I : iThesleff and Nieminenl . Il996l ) which, from 
a developmental biology perspective, greatly un- 
derestimate the role of the environmental com- 
ponent. Indeed, "development is phenotypic 
change in a responsive (plastic) phenotype due 
to inputs from the environment and the genome. 
[. . . ] Genomic and environmental factors have 
equivalent and potentially interchangeable devel- 
opmental effects, effects that depend as much 
on the structure of the responding phenotype 
as they do on the sp ecific inputs themselves" 
(|West-Eberhardl . l2003l ). 

The significant role of environmental factors 
in dete r minin g dental variation was reviewed by 
Butleil h98i ). who acknowledged that "dental 
phenotypic differences between populations do 
reflect genetic differences" (p. 288, Italics ours), 
but at the same time emphasized that the genes 
"that affect the dentition do so only through the 
mediation of a complex ontogeny in which the 
environment takes part" (p. 288). Different as- 
pects of the environmental influence on tooth 



development in human populations have been 
studied. One is the relationship between malnu- 
trition and tooth eruption, where malnourished 
childr en rev eal a delay in the eruption of their 
teeth (Boas, Il92 7: En wonwul. Il973l : iGarn et al.l 



19651 : iGaur and Kumarl . I2OI2I ) . Another way to 



study the environmental influence on dental de- 
velopment is to correlate tooth size and envi- 
ronmental stress. McKee and Lunz (|l990l ^ stud- 
ied the correlation between enamel hypoplasia 
occurrence, a marker of developmental disrup- 
tion, and dental size reduction: they found an 
association between dental reduction and the 
prese nce of enamel hypopla sia in the individ- 



ual. iGarn and Russell (|l97ll ) showed how ma- 



ternal health status can influence the tooth size 
of the deciduous and permanent dentition of the 
offspring. Other authors demonstrated that en- 
vironment plays a role ir i dete r mining too t h size 



(IKolakowski and Bailitl. Il98ll : iPerzigianl . Il984l : 



Townsend and Brown . 19791 ). Similar results 



come from the study of fluctuating asymmetry, 
which is considered a measure of stress and de- 
velopmental instability also influenced b y en- 
vironmental v a riables ( 
Parsons . 199d ). Kieser ( 



Carter and Hould . 



2011 



I992I ) , studying fluctuat- 



ing odontometric asymmetry, used alcohol con- 
sumption by the mother during the gestation pe- 
riod as an environmental variable, and found a 
high level of asymmetry in children of alcoholic 
mothers compared to childr en of non-al c oholic 
mothers. In another paper, (jKieser et al.l . Il997l l 



fluctuating asymmetry in children was consid- 
ered in relation to maternal obesity and smoking; 
again asymmetry was higher in children of obese 
and smoking mothers. The study of monozy- 
gotic twin pairs is another research approach 
to address how differences in the presence, size, 
and shape of teeth are linked to environmen- 
tal influence. "[...] Supernumerary tooth for- 
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mation is influenced not only by genetic fac- 
tors but als o by environmental a nd epigenetic 
influence s" (ITownsend et al.l |2012| . p. 6). For 
example, Townsend et al. ( 20051 ) found monozy- 
gotic twins to be discordant in the number and 
position of supernumerary teeth. Other stud- 
ies on monozygotic twins focused on the her- 
itability of discrete traits; in the first half of 
the last century, the prevailing view was that 
morphological traits on tooth crowns were un- 
der s trong hereditary control (jScott and Turner 
19971 . p. 161). For example, in the case of Cara- 
belli's trait, the most studied trait on upper mo- 
lar crowns, starting from the 1950s this view 
began to shift from a simple single-locus model 
to a more complex one, invoking multi-loci ge- 
netic control or environmental effects. As stud- 
ies on twins increased in number, it became ev- 
ident how trait expression often differs between 
identical twins (jScott and Turnerl 119971 . p. 162). 
Lastly, the effects of pathological conditions on 
dental development ( Atar and Korperich . 20 id ) 
have been described. Several pathological condi- 
tions {e.g. cystic fibrosis, AIDS, leukemia) cause 
a number of enamel and dentine defects, acting 
during tooth development, in the pre-natal and 
post-natal periods. 

All these studies discussed above have some 
limitations when considered from an evolution- 
ary developmental biology perspective. First of 
all, the majority of them concentrates on metri- 
cal traits; when discrete traits are considered, 
the focus is on characters of the whole den- 
tition, such as congenital absen c e and super- 



numerary teeth (ITownsend et al.l. l2005l) , or o n 



pathological defects ( Atar and Korperich , 20ld ) . 
In studies on environmental infiuence, normal 
individual variation which describes the mor- 
phology of a single tooth is usually not taken 
into account. Secondly, in evolutionary biology. 



the relevant variation is that occurring at the 
population level, since morphological differences 
among taxa often stem from evolutionar y pro- 
cesse s taking p l ace at the population level (jMavd . 
19821 : Ijernvali l200d ^. Few studies have taken 



into account environmental infiuences on teeth 
at the population level: among them, some were 
on pop ulations who m i grated into new environ' 



ments (Lasker, 19451 : Scott and Alexandersen . 



199^ ). but they failed to demonstrate plastic- 
ity in dental traits, probably because of th e 



kind of stressor tested ( Scott and Turner , 19971 ). 
Thirdly, even when some degree of plasticity is 
recognized {e.g. differences in trait expression 
between MZ twins), these studies lack a devel- 
opmental explanation of the mechanism through 
which this plasticity is produced. This diffi- 
culty is particularly evident in heritability stud- 
ies, where the issue is trying to distinguish be- 
tween genetic and environmental variance, as 
if one genotype were bound to always produce 
the same phenotype and, as such environmen- 
tal variance is a mere noise altering the cor- 
rect correspondence between genotype and phe- 
notype. In this framework, any attempt to ex- 
plain the developmental process becomes sterile. 
On the other hand, if we accept that variability 
on which evolution ary factors act is the prod- 



uct of development ( Salazar-Ciudad et al. . 20031 ) 



and that development is the result of the interac- 
tions between genetic and envir onmental inputs 
and the developing phenotype (jWest-Eberhardl . 
20031 ). the mechanism involved in the shaping of 



variability becomes central to our investigation. 

In recent years, thanks to remarkable advances 
in the field of molecular biology, it has been pos- 
sible to begin to understand the molecular ba- 
sis of this mechanism and how the developmen- 
tal process leads to the final phenotype. Stud- 
ies on experimental animals have provided in- 
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formation on genes expressed in different tis- 
sues during tooth development; tliis research, 
together with a broad scale analysis of morpho- 
logical variability in mammalian dentition, has 
led to the proposal of a specific model of devel- 
opment for multicuspidate teeth. The develop- 
mental unit f or cusp development is th e primary 
enamel knot ( Jernvall and June . 200d ). "a clus- 
ter of cells in the central part of dental epithe- 
lium, facing the denta. 1 mese nchyme" (jButleii . 



19561 : IVaahtokari et al.1 . Il996l l. The 



primary 

enamel knot acts like an embryonic signaling 
center and expresses the same genes as involved 
in the organogenesis of other ectoder mal organs 



such as feathers, nai l s, and glands (IBeil. 12009 



Jernvall and Theslefll . l2000l : I Jernvall and Jung , 
2000); the majority of these genes belong to 
one of the signaling pathways of four gene fam- 
ihes: BMP, FGF, SHH, and WNT. Further- 
more, the dimensions of the primary enamel 
knot determine the position and dimensions of 
the s e condary enamel knots (I Jernvall and Jung . 



20001 : iJernvall and Thesleffl . \200d ). from which 



the cusps originate. Lastly, secondary enamel 
knots express the same genes as expressed in 
the primary enamel knot, i.e., in the gene net- 
work there are no s pecifi c signals for each cusp 
( Jernv all and Jung!, bood). Starting from these 



clues, Jernvall (1995) developed a model for sec- 
ondary enamel knots initiation: a new enamel 
knot initiates when there are enough epithelial 
cells available to form them; the number of ep- 
ithelial cells is regulated by the expression of 
FGFs in extant enamel knots. The formation 
of cusp patterns is therefore a cascade of events 
beginning with the formation of the primary 
enamel knot: "the relative sizes, numbers, and 
heights of the cusps [in seal molars] are sym- 
metrically correlated, anterior and posterior, to 
a major central cusp, as if the development of 



these additional cusps depends on rate of growth 
relative to a threshold in time, beyon d which fur- 



ther developme nt does not o ccur" (jWeiss et al 
19981 . p. 390; Ijernvalll lioool ). The dimension 



of the primary enamel knot regulates the for- 
mation of subsequent secondary enamel knots 
and the dimensions of those regulate the initi- 
ation of other enamel knots (and cusps) later 
in the developmental cascade. Tooth morpho- 
genesis is therefore a self-o rganizing sy s tem o r 



19981 ). 



a Bateson- Turing Process (jWeiss et al.l . 
where "once a set of differentiation factors is 
established [. . . 1 fu rther exogenous s ignaling is 



no longer needed" ( Weiss et al. 19981 . p. 376) 



When the first enamel knot appears, the tooth 
primordium has all the information it needs to 
complete its development. This model for tooth 
development is usually referred to as the cascade 
mode for tooth morphogenesis, or the morpho- 
dynamic model. Today fairly accurate compu- 
tational models for tooth development are avail- 
able, taking into account both t he genetic net- 
work and the cellular interactions ( Osborn . [1993!, 



20081 : ISalazar-Ciudad and Jernvall I200I |2010|). 
Starting from the cascade model, these attempt 
to reproduce exta nt variability in tooth diver- 
sity. For example, ISalazar-Ciudad and Jernvall 
(|20ld ) developed a model accounting for the 
morphological diversity in seal molars. Follow- 
ing the suggestions of the cascade model, this di- 
versity is obtained by acting on very few initial 
inputs representing genetic information: one in- 
hibitor and one activator for the differentiation 
of enamel knots, and one parameter represent- 
ing growth factors, upregulated in enamel knots, 
regulating cell proliferation or differentiation. 

With this background, the aim of this paper is 
to test whether environmental inputs can affect 
dental morphology acting on the developmental 
cascade. Scott and Turner (I1997I ) suggested that 
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stressors such as undernutrition and infectious 
diseases, already known to disrupt dental devel- 
opment leaving marks on teeth, could affect den- 
tal morphology as well. The occurrence of den- 
tal enamel hypoplasia is a good marker for these 
kind of stressors, since it is known to be linked 
to physiological stress, such as malnutrition 
or systemic diseases, duri ng tooth de velopment 
(iGuatelli-Steinberg and Luk acs. 1999: lTen-Catel . 
19981 : I Skinner and Goodmanl . 1 19921 ^. The re- 
sponse of the developing tooth to such stress is 
a slowing down of its development and the visi- 
ble result on the tooth su r face i s an enamel hy- 
poplastic defect. Jernvall ( 200d ) suggested that 
studies on variability at the population level can 
be useful in testing hypotheses in development 
biology: our hypothesis is that, in a homoge- 
neous population, individuals with severe hy- 
poplasia, that is to say, individuals whose devel- 
opment was probably disrupted by an environ- 
mental stress, are more variable in their molar 
morphology than individuals with no hypopla- 
sia. The expectation we want to test is that a de- 
velopmental disruption, due to an environmen- 
tal stressor, occurring in a population increases 
variability in molar morphology. In a single or- 
ganism, a disturbance during the morphogenetic 
cascade should cause subtle "errors" in the mor- 
phological output. These changes will be more 
evident as the "errors" accumulate in the devel- 
opmental cascade, that is to say, in later devel- 
oping cusps. If we look at the whole population, 
where disruption has occurred at different times 
and with different intensity in each individual, 
the logical expectation is an increase in morpho- 
logical variability, particularly evident in later 
developing cusps. 



Materials and Methods 

The specimens analyzed derive from the "Fabio 
Frassetto collection" housed in the Museum of 
Anthropology at the University of Bologna. The 
largest part of this collection is made up of over 
600 skeletons recovered from the cemeteries of 
the town of Sassari (Sardinia, Italy) in the early 
20th century. For the vast majority of the skele- 
tons, individual information is available (name, 
sex, date of birth, date of death, job, and cause of 
death). This sample can thus be considered very 
homogenous in terms of genetic structure of the 
population and physical environment. Selection 
criteria required the presence of the maxillary 
dentition, with at least two molar teeth from the 
same side. Specimens where attrition precluded 
analysis of cusp morphology were excluded. The 
occurrence of enamel hypoplasia was employed 
here as marker of environmental stress {sensu 
lato). Enamel hypoplasia is known to be linked 
to physiological stress, such as malnutrition or 
systemic diseases, during derital development 
j Guatelli-Steinberg and Lukacs^. 'l99j|; Ten-Catj . 



1998 



Skinner and G oodman. 1992,). The pres- 



ence of li near enamel hypoplasia s was recorded 
following iBuikstra and Ubelakerl (j 19941 ). taking 
into account the severity with a two-grade scale 
on the basis of the breadth and depth of the 
hypoplastic bands (fig. ']^. Two groups were 
then identified, with different stress levels: the 
first group included individuals with no signs of 
enamel hypoplasia in any of the present teeth; 
the second group included individuals with at 
least one tooth with one severe or three faint hy- 
poplastic events. Only individuals with a clear 
pattern of presence/absence of hypoplasia have 
been included. The final sample included 75 
individuals (34 females, 41 males). Group 1 
(no hypoplasia) included 30 individuals (13 fe- 
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Figure 1: Linear enamel hypoplasia on the maxillary 
dentition: on the basis of the breadth and depth of the 
bands, severe and faint grade s are distinguished (following 



bangs, severe and tamt grade s 
iBuikstra and Ubelakeilll994l ) 



males, 17 males); group 2 (severe or diffuse hy- 
poplasia) included 45 individuals (21 females, 24 
males). For the sake of convenience, from now 
on group 1 and group 2 will be referred to as 
"non-stressed" (NSG) and "stressed" (SG), re- 
spectively. For each individual of each group, the 
variability in the expression of maxillary molar 
cusps was observed. Expression of the relative 
development of the cusps was recorded in a semi- 
quantitative way us ing the ASU-DAS standard 



([Turner et al.l . ll99lh . Discrete traits considered 



were: metacone, hypocone, metaconule, Cara- 
belli's trait, and parastyle. Each trait was scored 
on a scale from to 5 (metacone, hypocone 
and metaconule), to 6 (parastyle), or to 7 
(Carabelli's trait), where is the absence of the 
trait. A two-sample Wilcoxon test was used 
to compare each discrete variable between the 
two groups. Hotelling test and two-sample 
MANOVA (Pillai test) were used for multivari- 
ate analysis. The p — values were corrected for 
repeated tests with Bonferroni's method. All the 
statistical a nalyses were carried out u s ing t he 
software R ( R Development Core Team . 20101 ). 



Results 

The analysis was carried out for each group 
(NSG and SG) with sex pooled; any further sub- 
division would have reduced sample sizes and af- 
fected reliability of the statistical analysis. Ta- 
ble 1 shows, for each molar and for each mor- 
phological trait (cusp), the comparison between 
NSG and SG in the frequencies of expression of 
each category. A general trend appears in most 
cases: in the SG, the variability of expression of 
the different grades (the frequency of each cate- 
gory) is greater than in the NSG and the higher 
grades of expression are more frequent. That is, 
in 10 of the 15 cusps considered (metacone on 

and M^, hypocone and metaconule on 
and M^, Carabelli's trait and parastyle on 
and M^), the mode of the distribution is con- 
stant, while the mean increases. This condition 
is suggestive of a directional increase in cusp ex- 
pression variability in the SG. These differences 
in the distributions are significant (two-samples 
Wilcoxon test, p— value < 0.05) for the following 
cusps: Carabelli's trait on each of the three mo- 
lars; metacone on the first, hypocone on the sec- 
ond and metaconule on the third molar. When 
the p — values are corrected using Bonferroni's 
method, the values under the threshold of signif- 
icance are those of hypocone on the second molar 
and Carabelli's trait on the first and second mo- 
lars. The variability for parastyle is too low to 
make sense of it using any statistical test, but 
it is worth noting that the only two parastyles 
observed in the sample belong to two individuals 
of the SG. As a further step in the study it was 
deemed suitable to perform a multivariate anal- 
ysis in order to take into account at the same 
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Table 1: p — values for the Wilcoxon test applied to compare the 
distributions of cusp development between the non-stressed (NSG) and 
stressed (SG) groups. The threshold of significance was set at 0.05; 
the values that maintain their significance after Bonferroni's correction 
(p — value < 0.0033) are in bold. 



M2 

Metacone (p = 0.037)*** (p = 0.431)ns (p = 0.289)ns 

NSG SG NSG SG NSG SG 



1 

2 1 1 

3 4 8 10 18 

4 12 9 35 39 15 16 

5 31 64 12 26 9 

W IP w 

Hypocone (p = 0.302)ns (p:^0.002)*** (p = 0.408)ns 

NSG SG NSG SG NSG SG 

16 5 11 16 

1 8 9 4 6 

2 3 11 5 3 

3 4 4 10 18 6 14 

4 15 22 14 17 1 4 

5 23 46 11 

w w w 

Metaconule (p = 0.362)ns (p = 0.580)ns (p = 0.008)*** 

NSG SG NSG SG NSG SG 

40 64 46 68 26 33 

1 5 3 4 5 

2 1 1 

3 1110 2 

4 1 

5 2 



Table 1: continued on the next page. . . 
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U 




1 





1 


5 





13 













1 


6 





1 
















7 





3 
































Parastyle 




na 


(p= 


0.416)ns 


{p = 


0.440)ns 




NSG 


SG 


NSG 




SG 


NSG 
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43 


73 


51 




73 


27 


41 
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5 






















6 























Riga et al. 2013 



8 



Environmental Influence on Dental Morphology 



time the grade of development of every cusp 
in each tooth or in each single individual and 
compare it between the two groups. This was 
done since the scope of the work is to address the 
overall variability in cusp expression of the mo- 
lar teeth, rather than the relative development of 
each single cusp. A principal component analysis 
considering, for each individual and each group, 
all three molar teeth and the grade of expression 
of each cusp (except for the parastyle because of 
its very low frequency in the sample) was carried 
out. Figure ^ shows scatterplots of the first three 
principal components in pairs, where individuals 
of the NSG and SG are indicated in different 
colors. The first three principal components ex- 
plain the 77% of the variability of the dataset 
{PCI = 39%, PC2 = 23%, PC3 = 15%). In 
Figure 2 the weight of each variable on principal 
components is also shown (black arrows). Vari- 
ation along PCI is linked to Carabelli's trait on 
M^, hypocone on and M^, and metaconule 
on M^; variation on PC2 is mostly explained 
by Carabelli's trait on and hypocone on 
and M^; lastly, metaconule and Carabelli's trait 
on explain most of the variation along PCS. 
The same analysis was also performed separately 
for each molar tooth. A similar pattern in the 
spatial distribution of the individuals belonging 
to the two groups emerges in both types of anal- 
ysis. Individuals from the NSG tend to cluster 
together in a relatively delimited area of the plot, 
whereas those from the SG are more broadly dis- 
tributed on the Cartesian plane, partially over- 
lapping the NSG. PCI values for the NSG vary 
between —3.19 and 0.42, while in the SG between 
-2.88 and 5.86; PC2 values between -1.62 and 
2.10 for the NSG and between -4.12 and 2.17 
for the SG; PCS values between —0.70 and 0.17 
(NSG) and -2.18 and 4.09 (SG). This indicates 
a higher variability in the SG as compared to 



Table 2: Multivariate comparison (Hotelling 
test) between SG and NSG, for the molars pooled 
(M^+M^+M^) and for each molar separately. Signifi- 
cance threshold is set to p = 0.05 (p = 0.0125 after Bon- 
ferroni's correction). 



Tooth 


p — value 




Hotelhng 




0.00059 


Ml 


0.00000 




0.00004 


M3 


0.00178 



the NSG, more easily detectable on PCI, which 
best distinguishes the two groups, followed by 
PCS. Also, individuals in the NSG tend to fall 
into the margin of the distribution, suggesting 
that the enhanced variability of the SG is not 
isotropic but directional. Once again, the di- 
rection toward which variability increases is the 
direction of more developed cusps (indicated by 
the direction of the arrows in Figure 2). In or- 
der to test whether the differences in the ob- 
served distributions were statistically significant, 
a Hotelling's test was employed; first with 
the three molars pooled, and then for each mo- 
lar separately. In the first case, parastyle and 
metaconule had to be excluded from the anal- 
ysis, because of their overall very low frequen- 
cies and the reduced number of cases due to the 
omission of strings with NAs data. On the other 
hand, in the analysis of each molar tooth (where 
reduced variables for each case made it possi- 
ble to omit fewer strings), it was possible to in- 
clude metaconule, but not parastyle, in the vari- 
ables considered. In both analyses, differences 
between NSG and SG are statistically significant 
(Tab. [2|), even after Bonferroni's correction for 
repeated tests. 
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Figure 2: Plots of the first three PCs for all the traits in upper molars (parastyle was excluded because of the 
low frequency in the sample). The arrows represent the weight of each variable on the PCs (a=metacone on M^; 
b=hypocone on M^; c=metaconule on M^; d=Carabelli's trait on M^; e= metacone on M^; f= hypocone on M^- 

f^; h= CarabeUi's trait on M^; i 
M^: 1= CarabeUi's trait on M^). 



metaconule on M^; h= CarabeUi's trait on M^; i= metacone on M^; j= hypocone on M^; k= metaconule on 
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Discussion 

The different degree of expression of some dental 
morphological traits in the two groups of a mod- 
ern human population (SG and NSG) seems to 
suggest the existence of a link between enamel 
hypoplasia occurrence and upper molar teeth 
morphology, in that both are the result of the 
same cause: a developmental disruption, likely 
due to a non-genetic disturbance such as mal- 
nutrition or a systemic disease (among the most 
common causes of enamel hypoplasia). In other 
words, an environmental disturbance affecting 
dental development has an influence on the fi- 
nal phenotype. Two aspects of the results are 
worthy of note: first, as we hypothesized, a de- 
velopmental disruption resulted in an increased 
variability in dental morphology; secondly, the 
increase in variability is directional, toward an 
expression of larger cusps. These two aspects 
will now be discussed in detail. 

Variability increase in the stressed 
group 

The first effect of the developmental disruption 
on the expression of upper molars non metric 
traits is an increase of variation in the SG. 
The increase of variability in the expression of 
morphological traits under a stressful condition 
has usually been attributed to two different 
factors. The first is related to the existence of 
a hidden genetic variability that is canalized 
by a buffering mechanism in normal develop- 
ment ( Waddington . 19591 ): if an environmental 
stressor exceeds a threshold, the buffering mech- 
anism breaks down and new genetic variability 
arises (iRutherfordl . hood ). However, there are 
few cases in which this buffering mechanism 
has been discovered (IRutherford and Lindquist . 



19981 : ISollars et all . I2OO3I ) and there are no 
clues indicating that a similar genetic buffer- 
ing mechanism may be active during dental 
development. The other possibility discussed 
in the literature is that new variability arises 
from the dev el opmen tal mechanism itself 



( Salazar-Ciudad . 200?! ). As mentioned 



m 

the introduction, the current developmental 
model proposed for the molars has, as a logical 
consequence, an increase of variability in a 
developmentally disrupted population. Indeed, 
molar tooth development is a morphodynamic 
process in which morphogenesis and induction 
mechanisms are inter d epend ent (jjernvalll . Il995l : 



Salazar-Ciudad et al. . 20031 ). In morphody- 



namic systems a small change in the inputs can 
be amplified during the morphodynamic process 
(developmental cascade) and have unexpected 
results on the final morphology, because the 
relationship between the parameters' change 
and morphology is not linear, but depends on 
the developmental network, on the phase of 
morphogenesis, and on all the other inputs 
involved in tooth development. Therefore, 
the effect we observed in our sample on the 
increase of morphological variability in the SG 
is compatible with this current view of tooth 
development. Development is an epigenetic phe- 
nomenon; therefore, if we wish to understand 
the mechanism responsible for our results, we 
must look at the epigenetic level, where dental 
development is re gulated by cellular and ge- 
netic interactions ( Salazar-Ciudad and Jernvall . 
H). Both these factors (cehular param- 



eters and gene expression) can be affected 
by environmental inputs. As for the first 
factor, environmental stress is known to pro- 
duce a downregulati o n of cellular metabolism 



( Hand and Hardewig . 19961 ). The effects of 
this downregulation can be multiple; for ex- 
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ample during hypometabolic states (such as 
starvation or lack of nutrients), downregulation 
can affect both energy production and energy 
consumption, and can inhibit macromolecule 
synthesis and turnover. In a similar situation 
it is likely that other parameters important 
for tooth development may also change, such 
as, for example, epithelial and mesenchymal 
growth and differ entiation rates. In this r egard , 



some authors ( Johnston and Gallant 



Katso et al.l. l200ll : Kozma and Thomasl . 




Schmelze and Halj . \200di ) have suggested that 



cell growth is affected by environmental and de- 
velopmental conditions and that, in some cases, 
it can cha i ige in response to nutrient levels. 



Kim et al.l (j2002l ) identified Raptor, a protein 
sensitive to nutrient levels participating in the 
regulation of cell size. As for the importance of 
gene expression, there is a vast literature on how 
environmental stressors can alter gene expres- 
sion, even just considering nutritional stressors. 
In bees, for example, queens and sterile workers 
differentiate from genetically identical larvae 
thanks to a different nutrition that proba- 
bly influences the expression of some genes 
throu gh DNA methylation (jKucharski et al.l . 
20081 ) . Also, in Drosophi la, gene expressio n 



is regulated by starvation ( Zinke et al. . 20021 ). 



In mammals, diet can influence the pattern of 
methyl ation and the stab le expr ession of some 



genes (jjaenisch and Birdl . l2003l ): nutrition is 
important during the gestational period and 
can affect the phenotype of the offspring, for 
example induc i ng cha nges in methylation in rats 



(jBurdge et al. 



20071 ). Most of the literature 
of the environment on gene 
on regulation by changes 
but it is not the only 



on the effects 
expression focuses 
in DNA methylation, 
possibility. Transposable elements can also reg- 
ulate gene expression under stressful situations 



( Slotkin and Martienssen . 200?! ). For example 
in mammals, apoptosis regulati on is participated 
by a LTR retrotransposon, (jRomanish et al. 
20071 ) and retrotransposon activation 
known to be stress-induced (Frucht et al 



IS 



1991 



199S 



Hampar et al.l. Il976l: iHohenadl et al 



Ruprecht and Lanzavecchial . l2006l : 



Sutkowski et al.l . |2004| ). In humans, under- 



nutrition can alter the growth trajectory of 
babies; also, their morphology and physiology is 
affected by the nutr iti onal state o f the mother 
teateson et alLbood ). ISomel et al.l (120081 ) found 
that differences in diet can produce differences 
in gene expression. They replicated human 
and chimpanzee diets in laboratory mice and, 
after two weeks, they analyzed the differences 
in gene expression in the liver and in the brain. 
They found that 4-8% of genes expressed in 
the liver changed their expression; considering 
that only 15% of the genes expressed in the 
liver are differently expressed between humans 
and chimpanzees, these data are impressive. 
Thus during dental development, from teeth 
primordia in intrauterine life up to calcification 
of the permanent teeth, the environment can 
act on all the factors listed above. Changes 
in one or more of those inputs can affect 
tooth development, producing a wide range of 
morphological variation. 

Directionality of variability increase 

The other, somehow unexpected, main finding of 
the analysis is related to the evident directional- 
ity in the increase of the variability. This is quite 
apparent in the PC plots, where individuals be- 
longing to the SG are scattered in the Cartesian 
plane not randomly with respect to the NSC in- 
dividuals. In particular, along PCI, NSC indi- 
viduals tend to fall mostly toward negative val- 
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ues of the first component, whereas SG individ- 
uals are uniformly distributed along the axis. 
Thus variation increases in the direction of more 
developed cusps (fig. 2), a pattern which is also 
the most common in the univariate comparisons 
of each cusp. This result seems to be in co ntrast 
with the results of McKee and Lunz ( 199d ). who 
reported a correlation between enamel hypopla- 
sia occurrence and dental size reduction. This is 
not the case, since the ASU-DAS system, used 
in this study, takes into account the relative de- 
velopment of each trait in the tooth and not the 
absolute dimensions. We then considered two 
options to evaluate the source of the directional- 
ity observed: (i) it is a mathematical artifact due 
to the presence of a left wall in the distribution 
of ASU-DAS grades; (ii) it is a true directional- 
ity, given by the particular stressor considered, 
which alters development in a specific direction. 

Directionality as a mathematical artifact 

In general terms, when comparing two distribu- 
tions with different mean values, it is reasonable 
to consider the distribution with a larger mean 
value as the effect of some directional factor. 
However, in some cases, it may happen that a 
non-directional increase in variability results in 
a directional change in the distrib ution, withou t 
any directional factor in action (Gould, 19961 ). 
This is the case for asymmetrical distributions 
with a left wall, a mode standing constantly near 
the wall and a right tail of variable length. The 
increase in the mean may be explained as an in- 
crease in variability which, because of the asym- 
metric distribution, is resolved at all in the right 
tail. In our sample (considering only statistically 
significant comparisons between NSG and SG), 
Carabelli's trait on and metaconule on 
seem to fall in this category. The mode of the 



•NSG' 



•SG" 



L 



Carabelli's trait on M2 



Carabelli's trait on M2 



Metaconule on M3 



I 



Metaconule on M3 



Figure 3: DifFerences between SG (left column) and 
NSG (right column) for Carabelli's trait on and meta- 
conule on A/''. The directional increase in the mean of 
the SG could be a mathematical artifact due to a non- 
directional increase in variability in an asymmetric distri- 
bution with a left wall (absence of the trait) and a mode 
standing next to the wall. 

NSG is (absence of the trait) and variations 
from this state are absent or very rare; in the SG 
the mode does not change and more individuals 
fall into the right tail, which is also longer than 
in the NSG (fig. 2^. The increase in variability 
is resolved at all in the right tail, because values 
lower than (the left wall) simply do not exist. 
Therefore, the directionality in changes in mor- 
phology in the SG might derive from some non- 
directional factor acting at the epigenetic level. 
This could mean that the environmental stressor 
considered acts on a wide range of epigenetic pa- 
rameters, without any specificity, randomly with 
respect to the effects on the developmental cas- 
cade. 
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Carabelli's trait on Ml 



■ ■■■■I. 



Carabelli's trait on M1 



1 2 3 4 5 

Hypocone on M2 



12 3 4 5 

Hypocone on M2 



Figure 4: Differences between SG (left column) and 
NSG (right column) for Carabelli's trait on and 
hypocone on A/^. In this case, the increase in the mean 
of the SG cannot be due to a mathematical artifact, since 
the mode changes also; the most probable explanation is 
the action of a directional factor on the developmental 
response. 

Directionality as real 

The likelihood of directionality mathemat- 
ical artifact, however, does not fit with other 
results. Again, considering statistically signifi- 
cant results only, Carabelli's trait on and 
hypocone on M'^ show a different pattern of dis- 
tribution in the comparison between NSG and 
SG (fig. 4). In these cusps the mean value in- 
creases in the SG but, unlike the former cases, 
the mode also increases in the SG. In particular 
for the Carabelli's trait on M^, the whole distri- 
bution in the SG appears to be shifted toward 
higher grades, with a decrease in the "grade 0" 
and an increase in all the other grades. 



This evidence seems to support an explana- 
tion based on directional factors acting at the 
epigenetic level. Since dental development is an 
auto-organizing system, with a few starting in- 
puts and a consequential cascade of events, with 
no apparent specific signals for different cusps, 
it is conceivable that the pattern of directional 
increase in the variability in the SG is produced 
by the same mechanism for all cusps. For this 
reason, we cannot rely on an explanation based 
both on a non-directional and a directional fac- 
tor. At the moment, a directional factor act- 
ing at the epigenetic level seems to us to be 
the most probable explanation, since it can ac- 
count for both types of distributions observed, 
whereas a non-directional factor fails to explain 
the distributions for Carabelli's trait on and 
hypocone on M^. If we accept this hypothesis, 
then we might expect that the stressor consid- 
ered acts not randomly on the variety of param- 
eters affecting dental development. Among these 
parameters, regulation of signaling during cusp 
development seems to be the most important 
source of variation among individuals, while cel- 
lular param eters best explain variation along th e 
tooth row ( Salazar-Ciudad and ■Ternvali hoid ). 
Considering only molecular signaling, cusps pat- 
tern is controlled by a game of interactions 
between signaling molecules regulating epithe- 
lial and mesenchymal cell proliferation, s i ich as 



FGF-4, BMP, and Sh h (I.Ternvall et al.l . Il994 : 



Vaahtokari et al. . 19961 ). These interactions reg- 



ulate the induction of new enamel knots, that 
is to say they produce the cusp pattern in a 
forming tooth. BMPs (Bone Morphogenetic 
Proteins), in particular, are associated with 
enamel knot apoptosis and its termination, for 
example indu cing inhibitors of c ell prolifera- 
tion like p21 



(Kassai et al, 



Jernyall et al. . 19981 ) and ectodin 



20051 ): studies on ectodin expres- 
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sion demonstrated that ectodin-deHcient mice 
have extra teeth and extra cus ps originated 



from a larger primary enamel knot ([Kassai et al 
20051 ). Therefore, an environmental stressor 



could produce the pattern of increased cusp size 
and number simply inhibiting the production of 
BMPs, or BMP-induced molecules. Studies on 
gene expression during dental development are 
needed to confirm or reject this hypothesis. 

As a final remark, we would like to con- 
sider the possible evolutionary significance of 
the kind of plasticity observed in dental de- 
velopment. It has been shown that the evo- 
lutionary history of mammalian dental pat- 
tern is characterized by repeated convergent 
or parallel evolution of new cu s ps such as the 
hypo c one (IHunter and Jernvalll . Il995l : iJernvalll . 



2000: Jernvall et al 



19961 ). whereas cusp loss 
is an uncommon event. That is to say that 
the developmental response to environmental 
stress appears to be in the same direction as 
the evolutionary trend in mammalian denti- 
tion. There are several possible explanations. 
For example, the developmental response could 
be simply the mark of a developmental con- 
straint biasing variation in a certain direction; 
a non-isotropic variation due to developmen- 
tal reasons may re sult in a directiona lity im- 
posed to evolutionlGouldl, Il977l . I2OO2I ). An- 



other intriguing possibility is that developmen- 
tal plasticity in dental development plays an ac- 
tive role in evolution. Developmental plastic- 



tion in evolutionary studies ( 


Aubret and Shine 


2009: Chaoman et al. 2000: 


3mlen et a 




2007: 


Gibson and Hognesd. 1996: 


Price et al. 




2003: 


Pried. 2006: Rutherford and Lindauist 




1998: 


Sollars et al.. 2003: Stauffer and van Snik Grav. 


2004: West-Eberhard. 200c 


i). For 


exam- 



pie, plasticity in certain morphological traits 



could allow survival in difficult times; phe- 
notypic changes, with the time, may ac- 
commodate and le ad to a genetic change in 
the same direc t ion (Pigliucci and Murren . 20031 : 
Pigliucci et al.l . l2006l ) . Scholars are still debating 
this t opic and there is no consensus (|de Jona . 
2005 : Pig liucci et al.l . I2OO6I ). but it is certainly 
an interesting avenue for future investigations — 
also with special reference to dental evolution — 
in particular because nutritional stress as a result 
of environmental pressures is common in many 
species and may have played a significant role in 
human and mammalian evolutionary history. 



Conclusions 

The results of our analysis seem to support the 
hypothesis that an environmental stress can lead 
to changes in dental morphology. In the group 
with high levels of enamel hypoplasia, an in- 
crease in cusp dimension variability was iden- 
tified. Further, the increase is not isotropic, 
but directional: individuals with hypoplasia tend 
to have relatively larger cusps than individu- 
als without hypoplasia. These results open up 
new possibilities for future research; in par- 
ticular it would be interesting to inquire into 
the molecular mechanisms acting at the epige- 
netic level producing the developmental plastic- 
ity in molar morphology, for example testing 
the environmental infiuence on the expression 
of genes for BMPs or BM P-induced inhibitors 
such as p21 and ec todin (see Jernvall et al. 19981 : 
Kassai et aDl2005l ). Our results confirm, once 
more, how dental development is a complex phe- 
nomenon: it depends on inputs from the gene 
network and from the environment acting on the 
developing phenotype. Considering teeth as a 
strictly genetically determined biological struc- 
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ture, can be misleading in certain kinds of stud- 
ies. 
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